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Recently, there have been hints for dark matter (DM) annihilation in the galactic center to one or more photon 
lines. In order to achieve the observed photon line flux, DM must have a relatively large effective coupling to 
photons, typically generated radiatively from large couplings to charged particles. When kinematically acces- 
sible, direct annihilation of DM to these charged particles is far too large to accommodate both the DM relic 
density and constraints from the observed flux of continuum photons from the galactic center, halo and dwarf 
galaxies. We discuss three exceptions to these obstacles, generating the observed line signal while providing the 
correct relic density and evading photon continuum constraints. The exceptions are (i) coannihilation, where the 
DM density is set by interactions with a heavier state that is not populated today, (ii) forbidden channels, where 
DM annihilates to heavier states that are kinematically blocked today, but open in the early Universe, and (iii) 
asymmetric DM, where the relic density is set by a primordial asymmetry. We build simple models to realize 
these scenarios. 



I. INTRODUCTION 

Dark Matter (DM) is one of the primary pieces of evidence for physics beyond the Standard Model (SM). Although its particle 
physics nature remains a mystery, in many theories DM is a weakly interacting massive particle (WIMP), whose its relic density 
is determined by annihilation to SM particles with a weak-scale cross section. The same annihilation processes that sets its 
density can give rise to observable photon signals that can be observed in the Universe today. Radiation from DM annihilation 
into electrically charged particles produces an additional component to the continuous 7 spectrum. More strikingly, DM can 
annihilate directly into 77, •yZ, or 7/1 through processes with charged particles in loops, producing monoenergetic 7 lines. Since 
7 lines are not easily mimicked by astrophysical backgrounds, they are a "smoking gun" signature for DM. 

Recently, several groups 1 1-4 J have reported a 7 line spectral feature at E 1 sa 130 GeV in publicly available data from the 
Fermi Large Area Telescope (LAT) [5|. It is an exciting possibility that WIMP DM may explain this signal. The required 
annihilation cross section to 77 is (cv) 77 « 10~ 27 cm 3 /s depending on the DM profile, an order of magnitude smaller than 
needed for the relic density 0] . Moreover, Ref . (4] |6] reported evidence for an additional 7 line at E 1 « 111 GeV; the pair 
of lines is kinematically consistent with DM with mass m x 130 GeV annihilating to both 77 and 7Z, or m x ss 140 GeV 
annihilating to both 7Z and 7/1, as in the model of GJ-[9]. The Fermi collaboration has not confirmed these results, and null 
results from their most recent 7 line search [ 10] are in tension with the cross section required to explain this signal. 1 Also, the 
77 cross section is consistent with null searches for 7 lines from dwarf galaxies ifTTI . Astrophysical sources have been suggested 
to explain the signal lfl^ll4l . though Refs. |4][T5][16] suggest that the signal may prefer a DM interpretation if it is not due to 
instrumental effects. In any case, further analysis with more data is required before claiming a definitive discovery of DM. 

Since the electric charge of DM is zero (or extremely tiny IfTTI ). 7 couplings to DM arise radiatively. For example, in many 
WIMP models, DM couples to SM charged particles {e.g., fermion pairs // or WW) through weak-scale mediators, giving a 
tree-level annihilation channel \X ~~ * If or WW. In this case, \\ — > 77 arises at one-loop, through virtual charged SM particles 
as shown in Fig. [T] From the effective field theory point of view, these one-loop processes lead to a dimension six operator 
FfiuF^ u for scalar DM and dimension seven operator XJsX^nvF 111 ' f° r fermionic DM, where F^ v is the electromagnetic 
field strength 1 18 -20 1. Therefore, one expects the enhancement of DM annihilation to charged SM states over 77 to scale as 

{(tv) f f.wwK av )ii ~ I 71 "/") 2 ~ 1q5 - (!) 

Some enhancement of the photon signal can be achieved by placing the charged virtual particles and DM in SU(2)l multiplets 
(see, e.g., Ref. [21 ]), though the ratio remains large. 

Due to the large ratio in Eq. ([TJ, a DM explanation for the 7 line signal faces two main obstacles. First, Fermi LAT observations 
place strong constraints on DM annihilation to // or WW from the 7 continuum, at the level of (av) ffww ~ O(few) x 
10 -25 cm 3 /s, depending on the final state particles iflOl l22l l23l . As a result, one naively expects an upper bound on the 77 
annihilation cross-section (cru) 77 < O(few) x 10 -30 cm 3 /s, which is well below what is needed to generate the Fermi signal. 
Indeed, the neutralino interpretation of the line signal has been disfavored by these arguments ll24l - |251 . 



Note the Fermi collaboration searched for lines using all sky gamma-ray maps 1 1 1 , and the results are obtained for |b| > 10° plus a 20° X 20° square at the 
galactic center, using Pass 6 processing. However, in the analysis by Ref. |T|, the search regions were optimized for DM signals and Pass 7 was used. 



FIG. 1: WIMP annihilation to charged SM final states (Left), e.g., fermions // or WW, generates annihilation to 77 at one-loop (Right). 



Second, the total annihilation cross section in the early Universe must be (<jv) ~ 3 x 10 -26 cm 3 /s to generate the observed 
relic density. For (ctv) 77 ~ 10~ 27 cm 3 /s, according to Eq. ([TJ, one expects XX ~ * ff or WW to be far too large, giving a 
relic density much smaller than observed. Even if tree-level annihilation is p-wave suppressed, the additional O(10) suppression 
from the DM relative velocity (squared) is not sufficient to avoid depleting the DM relic density. 

In addition to the dimension six or seven operators just discussed, fermionic DM may couple to photons through a dimension- 
five magnetic dipole operator XP^X^ixv or electric dipole operator X^^lbX^uu, where a^ v = — i[7 M , 7"]/2. This type 
of DM can be found in models where DM is a composite state [27-35], and was considered recently in connection with the 
Fermi line signal 11361 . Dipolar DM encounters similar challenges in explaining both the line signal and relic density, since the 
dipole operator mediates XX ff as we H as XX 77- F° r the magnetic dipole case, fixing (crw) 77 = 1CP 27 cm 3 /s gives 
(crv) f t > 10 -25 cm 3 /s, which gives a too-small DM relic density. In the electric dipole case, XX ~^ ff is p-wave suppressed, 
and the relic density is too large, unless there are additional annihilation channels. Furthermore, Dirac DM models with such 
large dipole interactions are excluded by direct detection experiments 11371 . 

So far, we have seen both the relic density constraint and the continuum photon bound strongly disfavor simple WIMP models 
for enhanced 7 line signals. To alleviate these tensions, we have to consider extensions to the simple WIMP models with 
designed features to enhance the 77 signal iTSTl [36l 13814521 . 

In this paper, we discuss three generic scenarios that are exceptions to these constraints, allowing for a large 77 annihilation 
rate while annihilation to fermions is suppressed compared to Eq. ([T}, both in the early Universe and in the galactic halo today. 
The three exceptions are: 



Coannihilation: The relic density is set by X1X2 ff, where xi is DM and xi is a next-to-lightest state nearby in mass. 
Annihilation to // is suppressed during freeze-out by the Xi~%2 mass gap, giving the correct relic density for 0(10 GeV) 
splitting. No annihilation to // occurs today since xi decays to xi an d is not populated. 

Forbidden channels: DM annihilates to charged fermions FF that are slightly heavier than the DM particles themselves. 
Due to the high velocity tail of the DM distribution, annihilation occurs in the early Universe, setting the relic density, but 
is kinematically forbidden today. 

• Asymmetric DM (ADM): The relic density is set by a primordial DM asymmetry, where a large annihilation rate xx^ ff 
is quenched by the DM chemical potential. After freeze-out, the asymmetry is washed out by DM particle-antiparticle 
oscillations due to tiny DM number-violating mass terms, xx^ ~* 77 annihilation can occur today with a large rate, while 
XX* ~~ * ff can be p-wave or chirality-suppressed. 

In the remainder of this work, we study in detail several minimal DM models as examples to illustrate each of these mechanisms. 
In each case, we show that an enhanced 77 annihilation rate can naturally be reconciled with the observed relic density and 
present 7 continuum constraints. 

In Sec. [II] we discuss coannihilation, presenting two models: (i) magnetic dipolar DM, and (ii) coannihilation with charged 



partners, which generates DM coupling to 77 at dimension seven. In Sec. Ill we consider a model with forbidden channels, 
and we derive the mass gap between DM particles and charged states required for the correct thermal relic density. In Sec. IV 
we present a scalar ADM model and discuss the ingredients necessary for generating the 7 line while remaining consistent with 
other constraints. Our conclusions are summarized in Sec.|V] We focus in this paper on models needed to explain the 130 GeV 
line, though we emphasize that our results are easily generalized to the case of multiple lines. 



II. COANNIHILATION 



In coannihilation scenarios, DM freeze-out is dominated by annihilation with a next-to-lightest state that is nearby in mass. 
For concreteness, we consider X1X2 — > ff, where xi is the DM, xi is the nearby state, and / is a SM fermion. We assume that 
the X1X2 coannihilation channel is dominant in the early Universe, while direct X\X\ annihilation is suppressed. If the mass 
splitting Am = m 2 — mi is comparable to the freeze-out temperature TV, coannihilation can provide a natural framework for 
enhanced 7 signals from thermal DM: 
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FIG. 2: Magnetic dipolar DM xi annihilates to 77, 7Z, ZZ (Left), while // occurs by coannihilation only with \2 (Right). 



• In the early Universe, the thermally-averaged coannihilation cross section is suppressed by a Boltzmann factor 
exp(— Am/T). For Am ~ TV, the coannihilation rate becomes moderately suppressed, requiring larger couplings to 
reproduce the correct thermal relic density. 

• In the present Universe, \2 is not populated, and therefore X1X2 — > ff does not contribute to any indirect detection 
signals. However, direct annihilation XiXi ~~ * 77 can occur, and the rate can be enhanced due to the large couplings 
required for thermal freeze-out. 

Ultimately, within a given model, there will exist a preferred parameter region for Am and couplings that can simultaneously 
explain the relic DM density and the observed 7 signal. In this section, we first discuss some preliminaries for computing the 
DM relic density, closely following Ref. Il53l . and then we consider specific models in parts A and B. 

Similar to single species freeze-out, the relic DM abundance for a general coannihilation scenario is computed by solving a 
Boltzmann equation 



(a cf {V)(n 
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K q ) 2 ) 



(2) 



where n x = J7 n Xi is the total Xi density. In writing Eq. |2) in terms of only n x , we assume the individual densities n Xi are in 
chemical equilibrium due to rapid Xif Xjf an d Xi Xjff processes, such that 



<? 4 (l + A 4 ) 3/2 exp(-xA l ) _ 



ffcff 



(3) 



We have defined x = mi/T, Aj = (m, — mi)/mi, and g c ff = J7 + AJ 3 / 2 exp(— ccAj), with degrees of freedom for 
Xi- The thermally-averaged effective cross section is (a e gv) = J7 . nrj(aijv), where <Xy is XiXj annihilation cross section 
and its thermal average is 
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The DM relic density today is given by 



1.07 x 10 9 GcV~ 



1/2 
9* m Pl 



dx 
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(5) 



where mp] ps 1.22 x 10 19 GeV is the Planck mass and 5* is the number of degrees of freedom in the thermal bath during 
freeze-out. The freeze-out temperature Tf = mi/xf is obtained by solving Xf = In (0.038 ffefr m i^Pl (& e sv) /y/g*Xf), which 
can be done iteratively. Alternately, one can directly solve Eq. Q numerically; for the cases we consider below, we find that the 
agreement with Eq. |5]l is better than ~ 1- 3% depending on the mass splitting. 

Now, we discuss two models which give rise to the Fermi line signal and a correct relic density with the coannihilation effect 
in the early Universe. 2 



2 To be clear, our models rely on the mass splitting Am to suppress (o" e ff v), which is dominated by large X1X2 and X2X2 annihilation cross sections. This is 
distinct from models where X1X1 annihilation is itself too large, and (<j c sv) can be suppressed by l/g c s by having a "parasitic" species \2 that does not 
annihilate strongly (see, e.g., I54II55I ). 
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FIG. 3: The mass splitting Am = mi — m\ required for Q xi h 2 — 0.11 (Left) and annihilation cross sections for 77 (solid),7Z (dashed) 
and ZZ (dotted) (Right) with respect to the dark matter dipole magnitude fiB- We take mi = 130 GeV and the nuclear magnetic magneton 
UN ~ 0.161 GeV" 1 . The vertical (horizontal) band on the left (right) panel indicates (cru) = (1.27 ± 0.32) x 10~ 27 cm 3 /s JT]. 



A. Magnetic dipolar dark matter 

Although the electric charge of DM must be zero or very small, DM can possess a sizable electromagnetic interaction through 
an electric or magnetic dipole moment E71 - I291 . As we show, magnetic dipolar DM can account for the Fermi 7 signal, and 
coannihilation plays an essential role in achieving the correct DM relic density. 3 We consider a Dirac fermion \ coupled to the 
hypercharge field strength through a magnetic dipole interaction, with Lagrangian 

jSf = ixH + m B xX + ^ (X C X + XX C ) + ^X^XB^ , (6) 

where \ c — ~ *7 2 X* is the charge-conjugated \ field. We have two mass terms: a Dirac mass mo, and a Majorana mass rriM, 
which splits \ mt0 two Majorana fermions xi,2 with masses = \mo ± ttim |. Taking m\ < 771,2, Xi is the DM. In terms of 
Xi,2> the magnetic dipole interaction becomes 

= ^~X2^XiF^ + |»^XiV, (7) 

where /j 7 = hbCw and nz = —/J-bsw, an d (fiw) is the (co)sine of the weak mixing angle. For Majorana states, only 
Xi X2 transition dipole moments are allowed. The photon and Z boson field strengths are and Z„„, respectively. 
DM can annihilate to 77, 7Z, and ZZ final states, through /j-channel X2 exchange, shown in Fig. [2] The cross sections are 



4 4 2 

jiZm{m2 

cr(xiXi -> 77> = i 2 1 2T2' ( 8a ) 



<r(XiXi -> 7^)« = ,00^4/0^2 , 0^2 ^2 < 8b ) 



1287rmf (2m 2 + 2m 2 - m|) 2 
Mz( m i ~ m|) 3 / 2 (2mim 2 + m?j) 2 
47rmi(m 2 + to 2 , — m fj) 2 



<HXlXl -> ^Z)7J - t _ 2 , ™2 ^2T2 ' ^ C > 



where mz is the Z boson mass. To explain the Fermi signal, we fix m\ — 130 GeV and (cru) XlXl _> yl w 1CP 27 cm 3 /s. For 
7712 ~ mi, the 7Z and ZZ cross sections are comparable; in particular, X1X1 ~^ 1% generates a second 7 line at an energy 
Ey = mi — mr z j{Am\) ~ 114 GeV, which may be indicated in the data [4|. We estimate the size of the f is required to the line 



3 The case of purely electric dipolar DM cannot explain the Fermi 7 line, since the coannihilation process X1X2 - ► // setting the DM relic density is p-wave 
suppressed [29]. Fixing the electric dipole moment to require ct(xiX1 ~~ ► 77)f ~ 10 — 27 cm 3 /s, the DM relic density is too large (even if Am = 0) unless 
additional annihilation channels are present. 
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FIG. 4: Coannihilation X1X2 ff' (Left), where xi > s DM, the coannihilating state X2 and mediator <j> carry electric charge, and /, /' are 
SM fermions. X1X1 — * 77 arises at one-loop (Right). 
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where /Xiv w 0.161 GeV -1 is the nuclear magneton. 

In the early Universe, coannihilation X1X2 — > ff provides the dominant annihilation channel, shown in Fig. [2] The cross 
section is 

<? (X1X2 ^ f f)v = ac } ^ 7 Q/ + — — 7 _ — z^T 1 „2 .2 r. — ITTv? — I ( 10 ) 



c w sw{s - m|) 



2c 2 H ,s 2 H ,(s-m|) 2 



where Q/ is the electric charge in units of |e|, T3J is the weak isospin, and ct is a color factor for fermion / (3 for quarks, 1 
for leptons). The X1X2 W + W~ cross section is 0(1%) of the total // cross section, and can be neglected. In addition, 
subleading X2X2 77) Z^, ZZ also impact the relic density, and the cross sections are obtained by switching mi and mi in 
Eqs. (g). 



Taking Eq. ( 10 1, we estimate the annihilation cross section cr(xiX2 ~^ ff) v ~ 1-7 x 10 -25 cm 3 /s for mi = 130 GeV and 
Hb ~ 3.6 x 10~ 3 /iat as preferred by the Fermi line signal. Clearly, a dipole which is large enough to generate the observed 
77 line will give rise to too large an annihilation to // both for the relic density and for continuum constraints in the halo if 
DM is a Dirac fermion (corresponding to ttim = 0). This problem is easily solved in a model where the components of the 
Dirac fermion are split. In this case, annihilation to fermions proceeds only via X1X2 ff, an d the annihilation rate will be 
suppressed by a Boltzmann factor exp(— Am/I/) with Am = m-i — mi. Since Tt 6 GeV for mi = 130 GeV, we expect 
Am ~ O(10) GeV for the suppression mechanism to work. 4 

We calculate the relic density of xi numerically by using Eq. Q. In Fig.[5](Le/f), we show the mass splitting between xi and 
X2 required for the correct DM relic density as a function of the DM magnetic dipole moment fis (solid). We can see that the 
relic density constraint requires a larger mass splitting for a larger [is as expected. For \ib preferred by the Fermi line signal as 
indicated by the vertical green band, the required mass splitting is ^ 7 — 10 GeV. Note for a large /is, the annihilation cross 
section to 77 becomes large enough to set the relic density without the presence of X2 in the thermal bath as indicated by the 
steep rise of the curve for hb ^ 7.5 x 10~ 3 /is ■ In Fig.[3](/?/g/!f), we plot annihilation cross sections for 77 (solid), 7Z (dashed) 
and ZZ (dotted) with respect to hb for mi = 130 GeV. 

In this model, X2 decays to xi promptly in the early Universe and it is not populated now due to the mass splitting. Thus, the 
model evades the continuum photon constraint. Since the preferred Am is too large for signals in direct detection experiments, 
the most promising way to explore this model is through the Large Hadron Collider (LHC) |36 37l l56ll57ll . 



B. Coannihilation with charged partners 

Next, we present another coannihilation scenario in which the coannihilating state X2 carries electric charge. To be concrete, 
we consider the following Lagrangian: 



-Sfmt = X2(gs + gpih)xi<t> + fig's + 9'pis)f<t> + h.c. 



(ID 



4 In contrast, Ref. [36] focused on dipolar DM with Am ~ O(100 keV), which is sufficient to avoid continuum and direct detection constraints. Although 
Am is too small to obtain the correct relic density by coannihilation, they argue that the Fermi line might be reconciled with the DM abundance by having 
both electric and magnetic dipole moments, or through momentum-dependent dipole form factors. 
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FIG. 5: Solid contours show masses (m2,m^) and couplings gs.p for {av) 11 — 10 _27 cm 3 /s, with gp — O.lgs (Left) and gp — gs 
(Right) and Q X2 — 1. Dashed contours show Qa m h 2 = 0.11, for different values of g = x/lg^l 2 + 3p| 2 - Gray region is excluded by 

n dul h 2 < 0.11. 



where /, /' are SM fermions, </> is a complex scalar, and gs.p, g'g p are scalar (S) and pseudoscalar (P) couplings. We assume 
X2 and carry electric charge Q X2 |e| = Q^|e| and are 5'C/(3)c-singlets. We take \i to be Majorana. 

The DM relic density is set by coannihilation X1X2 — > ff and X1X2 /'/> shown in Fig. [4] X2X2 - ► 77 j // aiso occurs 
through gauge interactions. The cross sections are 



V12V = cr{xiX2 -> ff)v = cr(xiX2 -> f'f)v 



|gp| 2 (|gsl 2 + |gp| 2 )(^i + ^ 2 ) ; 

87r((s _ m 2 )2+m 2 r 2 ) 



P22V = a( X 2X2 -> SM)v = a( X 2X2 -> 77> + ^ *(X2Xa -> /j> = (Q* 2 + (20/3)Q 2 2 ) ^ 

/ 



(12) 
(13) 



where the factor of J2 f N lQ) = 20/3 arises from the sum over all charged SM fermions except t, which is kinematically 
blocked. The partial widths entering are 



r(^//') = ^(bsl 2 + l<7pr 
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(14) 
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We work to lowest order in v, with the exception that we take s = (mi + m 2 ) 2 (l + v 2 /4) in Eq. ( |T2] > to properly account for a 
possible resonant enhancement ll53lk near the resonance, (012W) must be computed numerically according to Eq. Q. We neglect 
contributions to cr 2 2 from Z-exchange and WW final states, which depend on the specific SU(2)l x U (l)y quantum numbers 
of X2 and <f>. Finally, the total effective cross section is 



(CTcffw) = 2r!r 2 (a 12 v) + — ((J 22 v) 



(16) 



where r 12 are defined as in Eq. Q with gi = 2, g 2 = 4. The relic density is given by Eq. 

DM can annihilate directly into 77 at one-loop, shown in Fig. |4] generating the 7 line signal. The cross section, given in 
Ref. I58H601 . is 



(crw) 77 = a(xiXl -t ll)v 



a Qx2 m i 

647r 3 mj, 



'F + \g s \ 2 +F_\g P \ 2 y 



(17) 



We have defined 



F+ = 



a ± vab T . 1 _ . 

T -^/ 1 (a,6) + — I 2 {a,b) + 



2b±2Vab 
l + a-b 



I 3 (a,b) 



(18) 
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FIG. 6: Tree-level dark matter annihilation to heavy fermions in the forbidden case (Left). \\ — > 77 at one-loop (Right). 

where a = m\/m\, b = mj/m|, and the functions I n (a, b) are defined in l[58ll . In the m^ ^> m 12 Am = m 2 — m i limit, 
we have F + w (2 — it 2 ) and F_ w 2; however, for m^ ~ mi, 2, these approximations overestimate the 77 rate and we use the 
exact expression in our analysis. Also, we expect the rates for X1X1 ~* %1 t0 t> e comparable, although the exact prediction 
depends on the SU (2) 1 x U(l)y quantum numbers of \2 and <f>. 
In Fig. [5] we present numerical results for this model. 

• The solid curves show mass contours for (crv) 11 — 10~ 27 cm 3 /s, for fixed mi = 130 GeV and for different couplings 
gs, with gp = 0.1 gs (left panel) and gp = gs (right panel). The 7 line signal requires gs > 0(1) and m.2, m^ > mi. 

• The dashed contours show parameters giving the DM relic density Q,d m h 2 = 0.11, for different values of the SM fermion 
coupling g' = y/\g's\ 2 + \g'p\ 2 , with mi = 130 GeV and gs,p fixed by (cru) 77 . There is a clear resonance for w 
mi + m,2, with smaller values of </ and larger Am allowed. (The width is computed as a function of the given 
parameters.) 

• The gray region is excluded by J],j m /i 2 < 0.11. For Am < 5 GeV, X2X2 annihilation is not sufficiently Boltzmann 
suppressed, depleting \i provided xi an d X2 are in chemical equilibrium. (This holds for g$ p ~ 1, g' 3> 10" 7 .) 

Taking m 2 rj 135 GeV (corresponding to the edge of the gray region) gives ft^h 2 = 0.11 in a large region of parameter space 
(10 -7 <C g' 10 _1 , off-resonance) with little dependence on the other new physics parameters, since the relic density is set 
through electromagnetic interactions. That is, the new physics particles need not have large couplings to SM states, aside from 
their electromagnetic couplings. In any case, this coannihilation model presents a viable framework for explaining the DM relic 
density with an enhanced 7 line signal. 

III. FORBIDDEN CHANNELS 

The second exception occurs when all the virtual charged particles generating the DM coupling to photons have a slightly 
larger mass than the DM. Although the coupling between DM and the charged particles has to be strong to overcome the 
loop-suppression factor, the annihilation cross section to charged particles at tree-level is suppressed kinematically. During 
freeze-out, DM is non-relativistic and its typical velocity is ~ 0.3 c. If the charged particles have masses not far from the DM 
mass, annihilation to the charged particles can still proceed in the early Universe, albeit less efficiently. As a result, one is able to 
obtain the correct relic density despite the large couplings needed to generate a photon line. On the other hand, DM has a typical 
velocity ~ 10~ 3 c in the halo today so that the direct annihilation to the charged particles is kinematically forbidden, evading 
constraints from continuum photons. In Ref. |8|, this mechanism was used to generate enhanced DM annihilation to •yZ and 
7/1, with the forbidden particle as the t quark. Here, we investigate a different model with enhanced annihilation to 77, and we 
compute the required mass splitting between the forbidden states and DM to obtain the correct relic density and the Fermi line 
signal simultaneously. 

We proceed to estimate the relic density through annihilation to the charged particle pairs, xx ~> FF> where we use F 
to denote charged fermions heavier than DM. We begin by reviewing the discussion of ll53l . Since the velocity of the final- 
state particles is small, it is convenient to write the annihilation cross section in the form (crv) = (a + 6w 2 )u 2 , where v is the 
relative velocity of the initial-state particles, is the velocity of the final-state particles in the center of mass frame, and a and b 
characterize the s-wave and p-wave contributions to the annihilation cross section respectively as usual. 5 Note w 2 must present in 
the annihilation cross section because it is from the phase space of the final-state particles. Energy and momentum conservation 



5 The reader should not be confused with the mass ratios a, b defined in Sec. II Here, a, b refer to s- and p-wave cross sections only. 
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The important step in computing the relic density for the forbidden case is to evaluate the thermally-averaged annihilation cross 
section, given by 



(av) = ((a + bv 2 )v 2 ) = 



x 3 ' 2 
2tt 1 /2 



v 2 (a + bv 2 )v 2 e- v2x/4 dv, 



(20) 



2/'- 



where fj,_ = (1 — m^/mp) 1 ^ 2 . Note 2/i_ is the minimal velocity to activate the annihilation. The integral of Eq. (20l can be 
performed numerically. In the case of an s-wave cross section off resonance, an analytical result is possible 



(av 2 ) = 



where z = mp / m x and K\ is the modified Bessel function [53 1. The relic DM density of \ is 



1.07 x 10 9 GeV" 1 



1/2 

g* m Pl 



rOO 

Jx f 



' dx 



(21) 



(22) 



where as usual the freeze-out temperature is Xf = In (0.038gm, x mpi (av) / -Jg^xjj. 

Having reviewed the relic density calculation in the forbidden case, we consider a concrete example. We assume that the DM 
X is a Majorana fermion and it couples to charged fermions through a pseudoscalar mediator. Pseudoscalar couplings are needed 
for a scalar mediator case to obtain an s-wave annihilation to 77. The interaction Lagrangian is given by 



l 9 x 



0X75X + igF<t>F~f 5 F, 



(23) 



where F is a charged fermion. Since we need mp > 130 GeV, the only possible candidate for F among SM fermions is the 
top quark. If F carries SU (2) p quantum numbers, there are comparable annihilation cross sections to ZZ and Z7, while if 
F carries only hypercharge, the 77 channel will dominate as discussed for the dipolar DM case. In the limit mp > m x , the 
annihilation cross section to photons through an FF loop is 



(cry) 77 = a (xX -> 77)' 



1 a 2 g 2 x g 2 F Q\cpm 2 p 
4tt 3 (s - m 2 )2 + m^T 2 



arctan 



1 



^Jmp/m 



(24) 



where is the total decay width of <fi, Qp is the electric charge of F in units of |e| and Cp is its color quantum number. In this 
model, is a sum of the following decay widths 



mi lD7r * 1 



where the function Ay 2 (r) is given by Ay 2 (r) = 2r 1 /(r) with 

/(r) = 



arcsin y/r\ 
i/i„„ 1+^ 



Am 2 , 

r(0^77) = ^^ 

mi 2567r J 7 



for t < 1 
— in) for t > 1 



F „2 



g 2 p\A^ /2 (m 2 /4m 2 F )\ , 
(25) 

(26) 



In our numerical study, is computed as a function of the given parameters according to Eq. (25 1. 

Since \X ~^ 77 is dominated by the s-wave process, the thermally-averaged annihilation cross section (o"u) 77 equals (eru) 77 . 
The annihilation cross section to FF is 



(av)pp = a(xx -> FF)v = 



glgpCFm 2 



2n ( s _ TO 2 )2+TO 2 r 2 



f'2- 



(27) 



The relevant diagrams are shown in Fig. [6] 
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m (GeV) m<p (GeV) 

FIG. 7: Contours show the coupling g x gF (Left) and the heavy charged particle mass raj- (Right) required for (cru) 77 = 10~ 27 cm 3 /s and 
Q x h 2 = 0.11 as a function of the mediator mass in the forbidden case. We take m x = 130 GeV, Qf = cf = 1, g x — 0.5<7f (dotted), 
9x = 9f (solid), and g x = 2gF (dashed). All contours stop when g x gF ~ O(40). 



Using Eq. ( 24 1, we can estimate the magnitudes of coupling constants required to generate the photon line signal. When is 
far from 2m x , a large coupling constant g x gp ~ 0(4 — 10) is required for (crv) 77 = 10~ 27 cm 3 /s, depending on the mediator 
mass. While m^ « 2m x , the line signal can be enhanced dramatically. In this resonance limit, the annihilation cross section to 
photons can be approximated as 



(av) 



io-3w/ s f g ^fl {^m 2 f 260 ^ 2 fi^y. (28) 

7 l 2xl0- 4 jVl30GcV/ V m 4> J \ r / 



Therefore, the line signal can be enhanced significantly in the resonance case and the required coupling constants can be much 
less than 0(1). 

Next, we discuss the thermal relic density for \- Since \\ — > FF is dominated by the s-wave process, we only keep the a 
term in the expansion of (av) = (a + bv 2 )v 2 , which is given by 

1 g y F c F m x 



where s = 4to 2 /(l — i> 2 /4) with a minimal v as 2(1 — rn^/mp) 1 ^ 2 . In our numerical work, we take the thermal average on the 
whole annihilation cross section av 2 as in Eq. p0] ). This is important to calculate the relic density near resonance. We also have 
checked that one may take s — Am 2 F and use Eq. pT] ) directly if it is off resonance. 

To see how we can enhance 77 signals and obtain the DM density simultaneously in the forbidden case, it is suggestive to 
check the ratio of (<ti | ) 77 to (av) p p. Taking Qp = Cp = 1 and mp > m x , we have 

^^xio-^i^ri^jg. (30) 

(<Tv) F p v 2 (4m 2 - m 2 ) 2 + m 2 T 2 

We see that there are two effects can overcome the loop suppression factor and boost (crv) 77 with respect to (av) FF . The first 
is the phase space factor v 2 , for rap > m x , we have v 2 -C 1. The second boost factor is from a resonance effect. Since FF 
annihilation occurs at s w Amp, while 77 annihilation occurs at s m 4m 2 , the latter can be enhanced by a pole at ss 2m x . 
Both effects rely on forbidden channels. If mp <C m x , then u 2 ~ 1 and both FF and 77 annihilation have the same resonant 
enhancement because they have a same pole at m^ = 2m x . Therefore, a successful implementation of these enhancements 
relies on the mass gap between F and \- 

We present our numerical results for the forbidden case on two complementary panels of Fig. [7] In Fig.[7](Le/f), we show g x g F 
required for the DM relic density and (crv) 77 = 10 -27 cm 3 /s as a function of m^. For each point along the contour, the value 
of mp is given in Fig.[7](/?ig/if)- When m^ sa 2m x , \X ~^ 77 is enhanced and a small coupling constant is needed to generate 
the Fermi line signal. In this case, a relatively small mp is required to suppress annihilation to FF. It is interesting to note that 
mp has to be very close to 130 GeV to obtain the correct relic density when 77 is enhanced maximally. On the other hand, 
(av)pp is on resonance during freeze-out for ~ 2mp. Therefore, one needs larger mp to suppress the boosted annihilation. 
We also can see that the numerical result only has a mild dependence on the relative size of g x and gp. This is because the 
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FIG. 8: Scalar DM \ annihilating to SM fermions // (Left) and 77 (Right), where F is a new massive charged fermion. 



value of each coupling enters the calculation individually only through T^, which is only important near the resonance. The 
dependence is negligible for small m^, where the width is very narrow and it does not play a role. For > 300 — 350 GeV, 
the effect is more noticeable since more decay channels become kinematically accessible. In this model, the preferred value of 
m,F is ~ 130 — 165 GeV depending on parameters. With such heavy charged particles, it is clear that XX ~^ FF is forbidden 
kinematically in the galaxy today, and the model evades the continuum photon constraint. 



IV. ASYMMETRIC DARK MATTER 



Asymmetric DM (ADM) [61 1 provides a third exception for reconciling an enhanced 7 line signal with the observed relic 

x 



density. 6 We assume that DM x is a complex state carrying a U (l) x conserved charge, and that a nonzero \ chemical potential 



arises sometime before the freeze-out epoch, generating an asymmetry of \ over its antiparticle x 1 ■ In ADM freeze-out, the 
XX* annihilation cross section can be much larger than ~6x 10~ 26 cm 3 /s required for symmetric freeze-out. In this case, xx^ 
annihilation is quenched once x^ is depleted, and the relic x density is determined by the primordial asymmetry. This is similar 
in spirit to coannihilation, where the coannihilating state x* is suppressed by a chemical potential, rather than a mass splitting. 

DM annihilation can occur in the Universe today if the x asymmetry is washed out after freeze-out through x ^ X* oscilla- 
tions [68-71 ]. Particle-antiparticle oscillations are generic in a wide class of ADM models where, unless U(l) x descends from 
a gauge symmetry, one expects t/(l) x -breaking mass terms to arise, e.g., through Planck-suppressed operators. In this case, x 
and x^ are no longer mass eigenstates, and oscillations commence once the mass splitting between the real components of x is 
comparable to the Hubble expansion rate. 

We consider a model where x is a complex scalar with an interaction 

^nt = xF(g L P L + g R P R )f + h.c. , (31) 

where gL.R. are couplings, / is a SM fermion, and F is a new massive fermion carrying U(l) x with mass mp > m x . We 
assume /, F carry electric charge Qf\e\ — Qf\z\- DM directly annihilates to // at tree-level and to 77 at one-loop, shown in 
Fig. [8] Since one expects the former to be enhanced over the latter by 0(ir 2 /a 2 ), we must address how this model can generate 
the observed 7 line while avoiding 7 continuum constraints. 

The case of scalar DM provides a natural mechanism to suppress annihilation to //, thereby evading the 7 continuum con- 
straint. If x couples chirally, XX* ~> ff is p-wave or chirality-suppressed as a consequence of angular momentum conservation. 
Taking, e.g., gL = 0, we have 



487r(m x + mp)" 2 



°{XX" ff)v « * , (32) 



keeping only the leading terms in v 2 or m 2 . On the other hand, if g^ ~ gn 7^ 0, the leading contribution is s-wave and is not 
chirality-suppressed: 



\9L\ 2 \gR\ 2 m 2 F 



0"(XX T ~> ff)v ~ '7 2 A? 2 V (33) 



For example, in the case of / = r, the annihilation rate in the galactic halo today [v ~ 10 3 ) is 

10" 23 cm 3 /s x \g L \ 2 \gR? for g L ~ gR 



v v 1 6 x 10 ffl cm7s x \gR\ for g L — 



6 For early ADM works, see 1 62- 66|; for more recent works, see 1 67 1 and Refs. therein. 
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FIG. 9: Left: Matrix element A as a function of mass splitting uif — 
coupling \qr\ and mass splitting rriF — m x required for (crtj} 77 = 10 _ 
XX s T ^ cross sections for / = r case. Shaded region is excluded b} 





= 1(T 25 cm 3 /s 


/ / / r- 
' ' / 1 
' ' 1 ' 
1 '/ ' - 

' '/ ' 
/ J . 




= 1(T 26 cm 3 /s 


Jr / 1 J - 


i (O-V)yy = 


1(T 27 ci£/s__ ■ — 






- 1(T 27 cm 3 /s 


f Hdm excluded : 



0.001 0.01 0.1 1 10 100 1000 

m F -m x (GeV) 

m x , for m x = 130 GeV and m/ = 0. Right: Solid contour shows 
!7 cm 3 /s, for m x = 130 GeV and ql = 0; dashed blue contours show 
ADM relic density considerations (see text). 



taking mp ~ m x = 130 GeV. Clearly, 0(1) chiral couplings are consistent with 7 continuum constraints, while nonchiral 
couplings are much more strongly constrained. 
The cross section for ~ * 11 i s given by 



, ar) = a?Q } {\gL? + \g R \ 2 Y ]A? ^ 2 x 10 - 29cm3/s x Q 4 (bi|2 + \ gR f fW , i3 s s 



The matrix element A, computed in Ref. [7 1 for m / = 0, can be expressed as 

A = 2 - 2 log (1 - r) - 2T- 1 arcsin 2 (y^) , (36) 
where r = rn^/mp. The numerical value of A is shown in Fig. |9|(Le/f). Although A diverges logarithmically for r — > 1, we 



expect the analytical formula to break down when r w 1 — rnj /m^ since 0(rnj) terms have been neglected. 

In FigM(Right), we show numerical results for \X annihilation cross sections for m x = 130 GeV and g L = 0. The solid 
contour shows the coupling and mass splitting mp — m x required for (cru) 77 = 10~ 27 cm 3 /s to explain the Fermi 7 line 
signal. The required parameters are easily consistent with 7 continuum constraints on \X^ ~^ //• F° r example, taking f = r, 
the dashed blue contours show the XX* t t cross section, easily consistent with present constraints [10, 22, 23]. Note the 
cases with / = e, /i are even less constrained by continuum constraints due to the chirality suppression. 

Lastly, we discuss constraints from DM relic density considerations. ADM freeze-out in the early Universe requires a large 
annihilation cross section (crv) > 6 x 10 -26 cm 3 /s to deplete the symmetric x density, leaving behind the residual asymmetric 
component. Although xx^ ~^ If is suppressed today, annihilation is greatly enhanced in the early Universe in two ways: 
(i) the DM velocity during freeze-out is v ~ 0.3, enhancing the p-wave term in Eq. ( |32| ), and (ii) for mp — m x < 10 GeV, 
coannihilation becomes important. The total effective annihilation cross section is (see Sec.|IH) 

(<r cff v) = r 2 Hxx 1 ff)v) + 2r x r F (a( X F -»■ jf)v) + r 2 F (a(FF -»■ SM)v) (37) 

with coannihilation cross sections 

a[xF 7/> = *QU\9rf + \9n\>)rn x ^ ^ + %) 

om^(m x + m]f) mp 

where for FF annihilation into SM particles we include only the dominant electromagnetic terms, summing over 77 and all 
fermions except t. Considering the case where / = r, gL = 0, and m x — 130 GeV, the gray region in Fig. |9]is excluded by 
requiring (<r c gv) > 6 x 10 -26 at Xf — 25. That is, the 7 line signal is fully consistent with ADM freeze-out. Parameters where 
symmetric DM gives the correct relic density correspond to the border of the gray and white regions, and therefore DM must be 
asymmetric in this model to explain the 7 line signal. 

In addition, we require that x X* oscillations begin during or after the freeze-out epoch to wash out the DM asymmetry, 
giving rise to observable annihilation signals today. Therefore, the U (l) x -breaking mass splitting should be less than H(Tf) ~ 
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4 x 10 -8 eV; otherwise, DM is symmetric since the asymmetry is erased before freeze-out. If oscillations occur much later than 
freeze-out, the DM relic density today is fixed by an initial asymmetry of 0(3.5%) of the baryon asymmetry. On the other hand, 
if oscillations begin during or soon after freeze-out, residual annihilation occurs and larger DM asymmetries are required to give 
the observed tt dm 1170117X1 . 

V. CONCLUSIONS 

Recent analyses of Fermi LAT data have found evidence for a 7 line signal from the galactic center at E 1 ps 130 GeV, with 
potentially a second line around 111 GeV. If these signals originate from DM annihilation, the required annihilation cross section 
to 77 is (cv) 11 w 10 -27 cm 3 /s, relatively larger than in generic WIMP models. To explain an enhanced 77 rate, one requires 
large DM couplings to light charged states, e.g., fermion pairs // or WW, generating a 77 coupling at one-loop. One expects 
tree-level annihilation to // and WW to be enhanced over 77 by 0(ir 2 /a 2 ). Therefore, a WIMP interpretation of the Fermi 
line signal faces two obstacles: (i) annihilation to charged SM particles in the early Universe is too large to explain the DM relic 
density, leading to excessive DM depletion during freeze-out, and (ii) annihilation to charged SM particles in the galactic halo 
today is too large, in conflict with Fermi LAT constraints on the continuum 7 spectrum produced by final state emission. 

In this work, we have emphasized three exceptions to these obstacles. For each case, annihilation to SM particles in the early 
Universe and in the halo today is suppressed, allowing for large DM couplings and an enhanced 77 rate, while giving the correct 
relic density and satisfying 7 continuum constraints for DM mass m x sa 130 GeV. The three exceptions are: 

• Coannihilation: The relic density is set by coannihilation X1X2 — > //■ An 0(10 GeV) mass splitting between DM \i 
and the nearby state \i gives the right suppression to // to explain both the relic density and 77 rate. Annihilation to // 
is absent in the halo since \i is not populated today. One natural example is a DM transition magnetic dipole interaction. 
We also considered a simple model where DM coannihilates with a state carrying electric charge. 

• Forbidden channels: DM annihilates \X ~> FF< where F is a charged state slightly heavier than \- Annihilation to FF 
is kinematically forbidden in the halo today, but occurs in the early Universe due to the higher DM velocity. We obtain 
the correct relic density for mp ~ 150 GeV. 

• Asymmetric DM: Due to a primordial x asymmetry, DM annihilation \X ff becomes suppressed in the early Universe 
when the symmetric x^X* density is depleted, with the residual asymmetric component providing the correct relic density. 
If the asymmetry is later washed out (through oscillations), DM annihilation today can give an enhanced 77 rate, while 
// is p-wave or chirality-suppressed. 

We illustrated these exceptions using simple models, showing in each case that an enhanced 77 rate can be naturally reconciled 
with the correct DM relic density and 7 continuum constraints. Clearly, a broad range of model-building possibilities lies within 
the general framework of these exceptions, beyond the simple models we considered. 

Virtually all the models we discussed here have a common feature: the presence of new charged states with mass near the DM 
mass. Such charged states would be prime targets for LHC searches, and could play an important role in modification of Higgs 
couplings to 77. In addition to the line from 77 annihilation, there appears to be another, lower energy line around 111 GeV, 
which may be consistent with annihilation to jZ. Depending on the SU (2)l quantum numbers of the charged states generating 
the effective DM coupling to photons, this line may also arise from a similar process as the one that generates the 130 GeV line. 
We leave an exploration of these points for future work. 
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